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Theoretical  analyses  of  small-scale  bridging  of  crack  surfaces  by 
elastic-ideally  plastic  springs  are  presented  and  applied  to  the  study  of  the 
fracture  toughness  of  ceramics  reinforced  by  small  particles.  The  dependence 
of  toughening  on  particle  size,  concentration,  strength,  and  ductility  is 
explored,  and  relations  between  toughening  and  bridge  length  at  fracture  are 
given.  Available  experimental  information  is  examined  in  the  light  of  the 
analyses.  j  yytihb  (2a^. 

INTRODUCTION 


sn 


The  fracture  toughness  of  a  composite  consisting  of  a  ceramic  matrix 
containing  ductile  metal  particles  is  studied  theoretically  in  this  paper  on  the 
basis  of  the  hypothesis  (e.g.  Krstic  1983)  that  the  metallic  inclusions  toughen 
the  ceramic  by  the  mechanism  of  crack-bridging  (Fig.  1).  It  is  presumed 
that  the  faces  of  an  advancing  crack  in  the  ceramic  are  pinned  together  by 
intact  particles  for  some  distance  behind  the  crack  tip,  reducing  the  crack-tip 
stress  intensity  that  would  otherwise  occur.  The  length  of  the  bridged  zone 
and  the  amount  of  toughening  that  occurs  depend  on  the  breaking  strength 
of  the  particles,  as  well  as  on  the  deformations  they  may  undergo  before 
they  fail.  Clearly,  for  crack-bridging  to  be  effective  a  matrix  crack  should 
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tend  to  be  attracted  to  the  particles,  and  it  would  appear  that  a  necessary 
condition  for  this  to  occur  is  that  the  elastic  stiffness  of  the  particles  be  less 
than  that  of  the  matrix. 

In  a  preliminary  outline  of  some  of  the  present  results  (Budiansky 
1986)  a  bridging-spring  model  (Fig.  2)  for  the  partially  pinned  crack  was  de¬ 
scribed.  The  same  model  has  recently  received  extensive  mathematical 
study  from  Rose  (1986),  who  cites  a  variety  of  fracture  problems  to  which 
the  model  has  relevance.  In  this  paper  we  limit  our  attention  to  the  case  of 
small-scale  bridging,  in  which  bridge  length  is  small  relative  to  crack  length, 
specimen  dimensions,  and  distances  from  the  crack  to  the  specimen  bound¬ 
aries.  We  discuss  the  associated  bridging-spring  equations  and  their  solu¬ 
tions,  and  then  consider  the  implications  of  these  solutions  with  respect  to 
the  problem  of  particulate  reinforcement. 

BRIDGING-SPRING  MODEL 


Elastic  springs 

In  the  absence  of  springs  connecting  the  faces  of  the  plane-strain, 
Mode-I  crack  shown  in  Fig.  2,  the  stresses  near  the  crack  tip  would  have  the 
standard  distribution 

°ap  -  Kfap(0)/(27tr)1/2  (la) 

where  K  is  the  stress  intensity  factor.  In  accordance  with  the  assumption  of 
small-scale  bridging,  we  will  study  the  effects  of  the  springs  by  letting  the 
crack  be  infinitely  long,  and  assume  that  the  distribution  (la)  is  approached 
asymptotically  for  r  -»  <».  The  crack-tip  stress  distribution  will  then  become 

°ap  •  Kmfap(0)  Ob) 

for  r  -»  0.  We  intend  to  calculate  the  new  stress-intensity  factor  Km  and  the 
spring  forces  in  terms  of  the  "applied"  stress-intensity  factor  K,  the  bridge- 
length  L,  and  the  spring  characteristics.  For  the  case  of  linear  springs  we 
adopt  the  notation  of  Rose  (1986)  here,  and  write  the  spring  stress  as 


o  -  kEv/(l-v2) 


(2) 
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in  terms  of  the  crack-face  displacement  v,  Young's  modulus  E,  Poisson's  ratio 
v,  and  a  spring-stiffness  coefficient  k.  A  relation  connecting  Km,  K,  and  the 
stress  cr(L)  at  the  end  of  the  bridged  zone  follows  immediately  from  the  use 
of  the  Rice  (1968)  J-integral 

J  =  Jr(Wni  "  °opuo,inp)ds  =0  (3) 

taken  around  the  path  shown  in  Fig.  3.  (Here  W  is  the  strain-energy  density 
and  ua  is  the  displacement  vector;  the  horizontal  parts  of  r  lie  just  outside 

the  line  of  springs.)  The  result 

(1-v2)K2  d-v2)K2m  t  (l-vV(L) 

E  *  E  kE  1  ' 

(Rose  1986,  Budiansky  1986)  provides  an  equation  for  the  toughening  ratio 
XsK/Km,  in  the  following  sense.  Suppose  Km  represents  the  critical  stress- 
intensity  factor  for  crack  growth  in  the  matrix,  and  imagine  that  new  springs 
connecting  the  crack  faces  emanate  from  the  crack  tip  whenever  the  crack 
propagates.  Then,  if  we  set  the  peak  spring  stress  o(L)  equal  to  the  spring 
breaking-strength  oY>  crack  propagation  with  simultaneous  fracture  of  the 
last  spring  will  occur  for 


X 


1  + 


(5) 


The  J-integral  relation  (4)  has  an  energetic  interpretation.  During  the 
propagation  process  just  described,  energy  is  released  into  the  crack  tip  at  a 
rate  (per  unit  crack  advance)  given  by  the  first  term  on  the  right,  and  the 
second  term  is  the  energy  loss  rate  caused  by  fracture  of  the  last  spring.  The 
left  side,  therefore,  represents  the  energy  input  rate  provided  by  the  applied 
field. 


We  would  also  like  to  find  the  relation  between  X  and  the  bridge 
length  L,  and  for  this  purpose  we  will  have  to  determine  the  spring  stress 
distribution  o(x).  A  convenient  starting  point  is  to  write  the  integral  equa¬ 
tion 


(l-vVtx)  4(l-v2)KVx~  40-v2)  fL_ _ rvr+vr^ 


k£ 


E/5x 


xE 


Jo  o<x  •  )'°g( 


V|x  -x  '|  J 


dx  ' 


(6) 


for  o(x)  in  (0,L).  The  left  side  is  the  crack-face  displacement  given  by  (2)  in 
terms  of  the  spring  stress;  the  terms  on  the  right  are  the  displacement  due 
to  the  applied  K-field  (la),  and  the  reduction  of  this  displacement  due  to  the 
action  of  the  spring  forces  on  the  crack  faces  (Tada  et  al.,  1985).  The  substi- 


tutions 

4kx  ^  4kx '  o 

give 

•  x  * 

*  ’  * _  Kvnr 

f° 

g(s)  +  g(t)  log 
Jo 

fr,+rtY  - 

where 

vFti  r  -^s 

a  = 


4kL 


0) 


(8) 


(9) 


It  follows  from  (4)  that  the  ratio  X  «  K/Km  is  given  by 


X  =  [l-2g2(a)l",/2 


(10) 


The  numerical  solution  of  the  integral  equation  (8)  for  g(s)  is  described  in 
detail  in  Appendix  A.  The  results  for  X  given  by  Eq.(10)  are  shown  in  Table  I 
for  values  of  a  in  the  range  (0,30),  and  a  plot  of  X  versus  a  is  given  by  the 
top  solid  curve  in  Fig.  4. 

A  check  on  the  numerical  accuracy  of  the  results  was  provided  by 
another  formula  connecting  X  and  g(s)  that  follows  from  the  relation  (Tada  et 
al.,  1985) 

Km.K-N/|j’()^dx  (11, 


which  gives 
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X* 


(12) 


A  direct  proof  that  the  J-integral  formula  (10)  is  consistent  with  (12)  and  the 
integral  equation  (8)  is  given  in  Appendix  B.  The  results  for  X  given  by  the 
formulas  (10)  and  (12)  were  in  excellent  agreement  with  each  other,  and  also 
with  those  found  by  Rose  (1986)  by  a  different  procedure. 

Since  g(s)  »  Vs  for  small  a,  it  follows  that  X  -  (1  +a )  for  a  -» 0.  Also,  for 
a  -*  ,  Rose  has  demonstrated  the  asymptotic  relation 

X  -  xVa/2  (13) 

The  results  for  the  non-dimensional  spring  stress  g  ■  g(s;a)  are 
plotted  in  Fig.  5  against  the  non-dimensional  distance  (a-s)  from  the  "last" 
spring,  for  a  *  2,  5,  and  20.  Also  shown  for  comparison  are  the  results  for  g 
given  by  a  Wiener-Hopf  solution  for  the  limiting  case  a-oo  (Budiansky  1986, 
Rose  1986),  summarized  briefly  in  Appendix  C.  Note  that  g(s)  peaks  sharply 
at  s-a,  consistent  with  the  behavior 

g'(s)  -  log(o-s)  (14) 


for  s  a  that  follows  from  Eq.(8). 

The  model  will  now  be  extended  to  the  case  of  elastic-plastic  springs, 
in  order  to  provide  a  basis  for  the  study  of  particulate  toughening  when 
plastic  yielding  of  the  particles  is  important. 

Elastic-plastic  springs 

We  assume  now  that  for  increasing  v,  the  springs  obey  the  elastic- 
ideally  plastic  constitutive  relation 


o«kEv/(l-v2) 


for  v  £  Vy  ■ 


qy(i-v2) 

kE 


m  <Jy 


for  v  £  v y 


(15) 


f  r 
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Then  if  v(L)  at  the  end  of  the  bridged  zone  exceeds  vY,  the  J-integral  result 
(4)  generalizes  to 

(1-V2)K2  (l-V^K2,  (l-V^Oy  ^  v 

- — * — - +  — — —  +  2oY[v(L)  -  VY1  (16) 


If  failure  of  the  last  spring  is  now  presumed  to  occur  when  [v(L)  -  vY) 
attains  a  critical  plastic  value  vp,  then  the  toughening  ratio  becomes 


X  =  I  1  + 


kK2 1  VY 


If  plastic  yielding  (v  >  vY,  o  =  oY)  occurs  for  LY  <  x  <  L,  the  spring  stress 
distribution  in  (0,LY)  is  governed  by 


(l-v2)g(x)  _  4(l-v2)KV  x 
kE  =  EV3k 


4(1  -v 
xE 


cv  x  4(1 -v2)  ^  .  ,V1  ( -Tx+^nr  i  , 

- - a(x  )log  -- - - -  dx  ' 

k  nE  Jo  {y  |x - x T J 


In  non-dimensional  form  (see  (7))  this  is 


8(s)tJo  8(t)lo8(T^1)d,+8Y£'08(;S7?)d 


t  «  Vs" 


(0£s£aY)  (19) 


where 


aY  -  4kLY/x  ,  gY  -  g(aY) 


The  solution  of  (19)  for  prescribed  values  of  a  and  aY  is  described  in  Ap¬ 
pendix  A.  The  associated  values  of  vP/vY  follow  from  the  fact  that 
v(L)  ■  vP  +  vY  is  given  by  the  right-hand  side  of  (18)  evaluated  at  x«L.  Hence 


s 
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Vp/Vy  =  -1  +  — 

8y 


(21) 


and  the  corresponding  value  of  X  follows  from  (17)  as 

1-1/2 


‘T 

A  check  is  provided  by  the  formula 

,«Y 


(22) 


1 


"Jo  $*M*-*) 


-1 


(23) 


that  follows  from  (11). 

The  previously  discussed  results  for  elastic  springs  correspond  to 
<XY**a  and  vp/vy-0.  By  varying  ay  at  given  values  of  a,  solutions  were  found, 
by  interpolation,  for  a  number  of  prescribed  values  of  vp/vy.  The  results  for 
X  versus  a  are  shown  in  Fig.  4. 

Rigid-plastic  springs 

In  the  limit  of  very  large  vp ,  the  elastic  contribution  to  spring  energy 
may  be  ignored,  so  that  the  J -integral  result  (16)  becomes 


(1-v2)K2  (1-v2)K$, 


+  20yVp 


(24) 


and  the  toughening  ratio  is 


X  = 


1  + 


2EOyVp 

I&l-v2) 


1/2 


(25) 


The  result  for  X  in  terms  of  bridge  length  follows  from  (11)  with  o(x)«oy;thus 

.  -  OyVE” 

X=l+2-v/r— —  (26) 


fl  Oy^ 
v*  Km 
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PARTICULATE  TOUGHENING 


Particle-pinned -crack  model 

For  simplicity,  we  will  assume  spherical  particles,  and  suppose  that 

the  faces  of  an  advancing  plane  crack  are  pinned  by  the  particles  at  their 

equators  (Fig.  1).  Consistency  with  this  assumption  would  require  that  the 
matrix  area  on  the  bridged  zone  be  reduced  by  a  fraction  that  exceeds  the 
volume  concentration  c  of  particles.  However,  again  for  simplicity,  we  will 
suppose  that  the  matrix-area  concentration  of  each  crack  face  retains  the 
value  (1-  c)  that  corresponds  to  an  arbitrary  plane  cross-section  of  the  com¬ 
posite  material.  To  invoke  the  bridging-spring  model  of  Fig.  2,  we  will  iden¬ 
tify  the  spring  stress  c(x)  with  smeared-out  particle  stresses  cop(x),  where 
op(x)  is  the  average  particle  stress  at  x  —  that  is,  the  average  over  all  parti¬ 
cles  in  the  thickness  direction,  normal  to  the  plane  of  Fig.  1.  Next,  we  have  to 

get  an  estimate  for  the  effective  spring  constant  k. 

Effective  spring-constant 

Consider  first  two  elastic  half-spaces,  bridged  by  a  single  isolated 
elastic  particle,  that  are  pulled  apart  by  remote  tensile  forces  that  produce 
an  average  tensile  stress  op  in  its  equatorial  cross-section.  If  we  neglect  the 
difference  between  the  particle  and  matrix  moduli,  the  remote  displacement 
of  each  face  (Fig.  6a)  is  given  by  the  classical  smooth-punch  result 


v  = 


tta(l-Vm)gp 

2Em 


(27) 


where  Em  and  vm  are  the  Young's  modulus  and  Poisson's  ratio  of  the  matrix, 
and  a  is  the  particle  radius.  This  is  the  average  crack-face  displacement  for 
the  limit  c-»0.  For  arbitrary  c  we  will  write 


ave 


=  P(c)*a(1'v?)cr 

2Em 


(28) 


and  estimate  P  on  the  basis  of  the  circumferentially  cracked  cylinder  shown 
in  Fig.  6b,  with  outer  radius  b  chosen  to  satisfy  (a/b)2  =  c  in  order  to  meet 
the  assumed  crack-face  area  ratio.  (Here  vave  is  defined  as  the  displacement 
averaged  over  the  full  area  xb2  of  the  cylinder.)  For  this  configuration,  Tada 


ct  al.  (1985)  supply  the  following  approximate  relations  for  the  crack  stress- 
intensity  factor: 


K,=  o^”a(1  ~a/k  —  H(a/b) 

2(a/b)2 

H(a/b)  =  1  +  ^-(a/b)  +  |(a/b)2  -  .363(a/b)3  +  .73 1  (a/b)4 

2  o 


(29a) 

(29b) 


where  the  applied  stress  is  o  =  (a/b)2op.  The  energy  released  by  cutting  the 
crack  into  the  cylinder  in  the  presence  of  constant  a  is  given  by  the  integral 


b^O-v|) 

Em 


a'da' 


(30) 


and  this  may  be  equated  to  nb2ovave  .  The  resulting  expression  for  vavc  may 
then  be  compared  with  Eq.  (28)  to  yield  the  formula 


(31) 


Finally,  substituting  vave  into  the  defining  relation  (2)  for  k,  we  get 


Here  E  and  v  are  the  effective  elastic  constants  of  the  composite  material 
consisting  of  the  ceramic  matrix  containing  a  randomly  distributed 
concentration  c  of  particles. 

The  function  p(c)  defined  by  Eq.  (31)  is  shown  by  the  solid  curve  in 
Fig.  7.  The  dashed  curve  shows  the  approximation  provided  by  the  empirical 
formula 

P-(l-c)(l-Vc)  (33) 

Since  this  estimate  for  f)  is  explicit  and  simple,  we  will  use  it  henceforth. 


Toughening  ratio 

We  can  now  rewrite  the  J-integral  results  (4)  and  (5)  to  make  them 
applicable  to  a  ceramic  reinforced  by  elastic  particles.  The  left-hand  side  of 
(4),  which  comes  from  the  far  field,  stays  the  same,  but  now  E  and  v  refer  to 
the  effective  properties  of  the  composite.  The  first  term  on  the  right  came 
from  the  energy-release  rate  at  the  crack  edge,  and  so  we  should  change  E 
and  v  to  Em  and  vm  because  the  crack  edge  is  advancing  only  into  matrix 
material.  We  also  multiply  this  term  by  (1-c)  to  take  into  account  the  fact 
that  the  bridging  particles  reduce  the  length  of  the  advancing  crack  front  by 
the  factor  c.  This  gives 

K2(l-v2)  K£(l-Vm)(l~c)  (l-y2)o2(L) 

E  “  E*  kE 


Finally,  we  use  Eqs.  (32)  and  (33)  to  substitute  for  k  in  the  last  term,  and  set 
o(L)  *  oY  «=  cS,  where  S,  the  particle  strength,  is  the  magnitude  of  the  particle 
stress  Op  at  failure.  The  result  is 


K2(l-v2)  _  Km(l~Vm)(l~c)  |  n  S2ac(l-c)(l-V~ :)(l-v2) 
E  Em  2  Em 


(35) 


(Here  we  should  interpret  Km  as  the  root-mean-square  value  of  the  stress- 
intensity  factor  in  the  ceramic  along  the  crack  front.)  Then,  for  the  case  of 
particles  that  break  elastically,  we  get  the  modified  toughening  ratio  • 


where 


A  = 


K/K 


m 


V(0(  1-C) 


1  + 


n  S2ac(l-Vc”) 


1/2 


0)  = 


Ed-Vm) 

Em(l-v2) 


(36) 


(37) 


Similarly,  the  J-integral  results  (16)  and  (17)  for  elastic-ideally  plastic 
springs  can  be  adapted  to  apply  to  ceramics  containing  ductile  particles  that 
undergo  plastic  flow  before  they  fail.  For  the  configuration  of  Fig.  6  the 
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connection  between  average  particle  stress  op  and  vavs  up  to  failure  is  now 
presumed  to  be  like  that  shown  by  the  solid  curve  in  Fig.  8.  We  idealize  this 
function,  replacing  it  with  the  elastic-plastic  relation  given  by  the  dotted 
lines,  choosing  S  as  the  peak  particle  stress,  and  picking  vP  so  as  to  make  the 
areas  under  the  two  curves  the  same.  (Unfortunately,  little  information  is 
available  concerning  actual  op-  vave  relations,  so  the  idealization  is  only 
conceptual.)  Then  the  elastic-plastic  counterpart  of  (35)  is  found  by  multi¬ 
plying  its  last  term  by  the  factor  (l+2vp/vy),  and  the  modified  toughening 
ratio  becomes 


A_  K/Km 
Vco(l-c) 


1+ 


jj  S2ac(l-V~c) 


i+2^£ 

VY  j 


1/2 


(38) 


For  a  given  strength  S,  the  effect  of  pre-failure  ductility  is  clearly  substantial; 
however,  the  lowered  particle  constraint  associated  with  increasing  vP/vy 
may  be  expected  to  lower  S. 

Finally,  the  toughening  for  the  limiting  case  vp/vy  ->  «>  may  be 
calculated  by  modifying  the  J-integral  result  (24)  for  rigid-plastic  springs. 
Thus 


and  this  gives 


K2(l-v2) 

E 


K£(l-v2)(l^c) 

Em 


+  2cSvP 


(39) 


A_  */Km 

Vto(l-c) 


1  + 


2c  E^vp 
1-0  K£(l-v£) 


(40) 


Toughening  vs.  bridge  length; _ matrix  crack-growth  resistance 

With  X  replaced  by  the  modified  toughening  ratio  A  defined  above, 
and  the  use  of  Eqs.  (32)  and(33)  to  evaluate  k,  the  results  of  Fig.  4  for  tough¬ 
ening  versus  bridge  length  may  be  applied  to  particle-reinforced  ceramics. 
For  convenience,  the  curves  of  Fig.  4  are  rescaled  in  Fig.  9  to  show  A  versus 
(x2a/8),  which  now  equals 


Ps 


c  L 
(1— c)(l— Vc-* )  coa 


(41) 
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(The  symbols  in  Fig.  9  represent  experimental  data  to  be  discussed  later.) 
Besides  the  connection  these  curves  give  between  L/a  at  fracture  and  the 
toughening  ratio,  they  have  another  interpretation  as  resistance  curves  asso¬ 
ciated  with  matrix  crack  growth  that  is  not  accompanied  by  particle  fracture. 
Suppose  a  pre-existent  crack  in  the  composite  is  subjected  to  a  gradually  in¬ 
creasing  far-field  K  and  think  of  A  as  the  current  value  of  the  ordinate  of  Fig. 
9  during  this  loading  process.  At  A=1  the  crack  will  begin  to  advance  into  the 
matrix,  but  particles  will  remain  intact  as  long  as  A  is  below  its  critical  value 
for  composite  fracture.  If  we  now  let  L  denote  the  amount  of  bridged 
matrix  crack  growth  that  occurs  before  failure,  the  relation  between  K  and  L 
is  provided  by  the  curves  of  Fig.  9. 

For  vp/vy  ,  the  relation  between  toughening  and  bridge  length  is 
produced  by  changing  Eq.  (26)  to  read 

A=1+2\/i  IF  (42) 

Toughening  vs.  concentration:  bridge  length  vs.  concentration 

The  toughening  ratio  K/Km  generally  starts  out  as  an  increasing 
function  of  c,  but  its  rate  of  increase  may  drop  substantially  as  c  goes  up.  For 
the  case  of  very  ductile  particles,  this  is  illustrated  in  Fig.  10a  when  the  ratio 
of  matrix  and  particle  moduli  is  Em/Ep=5.  Here  Eq.(40)  has  been  used  to 
calculate  K/Km  versus  c,  for  several  values  of  Z«SVL0/Km,  where 


L0 


ZL  EmvP 

8  KJl-Vm) 


(43) 


The  parameter  co  defined  by  Eq.(37)  was  calculated  as  a  function  of  c  on  the 
basis  of  the  self-consistent  method  (Hill  1965;  Budiansky  1965). 

Similar  results  may  generally  be  expected  for  finite  vp/vy.  However, 
for  a  sufficiently  small  value  of  vp/vy,  K/Km  may  reach  a  maximum  at  a 
moderate  value  of  c.  (For  the  limiting  case  of  elastic  particles,  with  Em/Ep>l, 
Eq.(36)  implies  a  maximum  value  of  K/Km  at  c  <  .29) 

It  follows  from  Eqs.(40)  and  (42)  that  for  very  ductile  particles 
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L  4  (  0)  ^ 

Lo  (A+l)2l1-cJ 


(44) 


Since,  generally,  co<(l-c)  for  Eiu/Ep>l,  the  critical  bridge  length  is  a  decreasing 
function  of  c,  as  illustrated  in  Fig.  10b  for  vp/vy  -*  00  •  Furthermore,  it  can  be 
shown  that  L  is  a  decreasing  function  of  c  for  all  vp/vy. 


EXPERIMENTAL  INFORMATION 


Some  available  experimental  measurements  of  bridge  lengths  as  well 
as  toughening  in  several  composite  systems  will  be  examined  in  the  light  of 
the  present  analysis.  Krstic  et  al.  (1981)  tested  a  model  system  consisting  of 
glass  reinforced  by  aluminum  particles,  and  Sigl  et  al.  (1987)  report  results 
of  studies  of  AI2O3/AI  and  WC/Co.  The  toughness  K  and  the  bridging  zone 
size  L  were  measured  by  placing  cracks  into  the  composite,  imposing 
external  loads,  and  observing  the  stretched  metal  ligaments  between  the 
surfaces  at  final  fracture.  For  each  composite,  the  data  to  be  discussed  are 
for  a  20%  volume  concentration  of  particles. 

Table  II  lists  the  measured  values  of  K,  L,  and  the  average  particle 
radius  a;  bulk  elastic  constants  of  matrix  (Em,  vm)  and  particle  (Ep,  vp);  elastic 
constants  of  the  composites  (E,  v),  calculated  self-consistently;  matrix 
toughness  Km;  and  the  estimated  uniaxial  yield  strength  o0  of  bulk  particle 
material.  Based  on  these  data,  the  points  plotted  in  Fig.  9  show  the 
corresponding  values  of  A  and  p.  The  effective  values  of  vP/vy  that  may  be 
inferred  are  vp/vy  »  41  for  the  glass/ Al  system,  »  68  for  AI2O3/AI,  and  *  12 
for  WC/Co. 

These  inferred  values  of  vp/vy  may  be  used  in  conjunction  with  the 
toughening  formula  (38)  to  back-figure  effective  particle  strengths  S.  For  the 
particles  in  glass/ Al,  we  find  S  -  250  MPa,  which  gives  S/a0  -  3.6.  But  a 
much  higher  strength,  S  *  1500  MPa,  is  calculated  for  the  considerably 
smaller  aluminum  particles  in  the  AI2O3/AI  system.  This  gives  S/o0  ®  22,  and 
the  same  surprisingly  high  ratio  is  obtained  for  the  small  cobalt  particles  in 
WC/Co,  for  which  Eq.  (38)  gives  S  »  10  GPa. 

For  AI2O3/AI,  the  rigid-plastic  model  seems  appropriate,  and,  indeed, 
using  Eq.(42)  to  estimate  S  provides  a  result  only  slightly  lower  than  that 


-14- 


given  by  the  elastic-plastic  model.  Finally,  we  remark  that  the  use  of  this 
calculated  value  of  S  •  1500  MPa  in  Eq.(40)  to  compute  the  plastic  particle 
stretch  gives  vp  -  .3  pm.  Sigl  et  al.  (1987)  measured  stretches  -  1pm  in 
severely  necked-down  particles  before  fracture. 

CONCLUDING  REMARKS 

Our  theoretical  study  of  particulate  toughening  displays  the  effects  of 
particulate  strength,  deformation,  and  size  on  toughening.  In  the  case  of 
elastic  particles,  Eq.(36)  shows  that  for  a  given  particle  strength,  toughening 
is  an  increasing  function  of  particle  size.  The  same  is  true  for  plastically 
deforming  particles  (Eqs.(38,40))  if  we  assume  that  both  Vp  and  Vy  are 
proportional  to  particle  radius.  Large  plastic  deformation  of  the  bridging 
particles  before  fracture  can  have  a  substantial  toughening  effect,  and  high 
particle  strength  is  obviously  desirable.  But  the  actual  strength  and  ductility 
of  a  bridging  particle  are,  at  this  time,  difficult  parameters  to  assess,  since 
they  may  well  be  greatly  influenced  by  load  and  shape  asymmetries,  as  well 
as  by  interface  sliding,  all  of  which  tend  to  reduce  constraint  on  the  particle, 
promote  plastic  flow,  and  lower  its  strength. 

Confrontation  of  experimental  data  on  metal-reinforced  ceramics 
with  the  predictions  of  the  present  theory  seems  to  imply  suspiciously  high 
particle  strengths.  If  such  high  strengths  are  not  confirmed,  toughening 
mechanisms  in  addition  to  crack  bridging  may  be  operative  in  particulate- 
reinforced  ceramics. 
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APPENDIX  A 

Integral-Equation  Solutions 

Elastic  springs 

A  convenient  set-up  for  solving  the  integral  equation  (8)  is  obtained 
by  differentiation,  which  gives 


8  (S)  +  H  \  A  ^  dt  *  TT= 

2 Jo  V  *  2vs 


(Al) 


To  satisfy  the  original  integral  equation  (8),  we  have  to  enforce  it  at  just  one 
point,  say  a«s,  in  addition  to  satisfying  (Al)  for  s  in  (0,a).  Note  that  (8),  (12), 
and  (Al)  imply  that 

g(s)  -  for  s— *0  (A2) 


If  we  now  let 


then  (Al)  reduces  to 


a 

2sH '(s)-H(s)  +  sf  — dt=s 
Jo  t-s 


(A3) 


(A4) 


A  neat  way  to  proceed  is  to  write 


H(s)  -  A(s/a)  +  f(s) 


where  A-  Vag(a),  and  f(0)  -  f(a)  -  0.  Then 


A  (l + 1 }  +  o  l0‘  f2?  1  +  2sf'<5)  ‘ f<8)  +S1 TT  ' s 


(A6) 


and  with 


S  *  — (1-COS0) 
2 


the  Fourier-series  expansion 


f*  ]^iBsin(ne) 


(A8) 


»1 


it  appropriate.  Since 


cot  (nO) 


(A9) 


>■1 


Eq.(A8)  becomes 


-^na.tin  n8+^11^2*®2  ^aBcoi  n6  *  -(l-cos6) 


(A  10) 


■■l 


If  we  now  multiply  (A10)  by  sin(m6)  and  integrate  over  (0,a)  we  get 


CmA 


"  o  *« 


*  <*F.m  (m=l,2,...) 


(All) 


0*1 


where 


Cm*  (I-kxJEoj  +  ^  TB 


Fm  *  — 


1_ 

m 


(m  odd) 


m 


m2-l 


(m  even) 
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m 


2  2 

m  -  n 


(m+n  odd) 


m(m2  -  n2  -1) 


(m2-nl  -l)a  -4n2 


(in -hi  even) 


tad 


TB.£(l-co.e)J  log(i^2|2],,.(m8)de 


was  evaluated  numerically.  The  substitution  of  (A5)  into  the  undifferentiat¬ 
ed  integral  equation  (8),  and  its  assertion  at  the  point  s-a  leads  to 


PnA 


(A12) 


where 


Pn  *  —  + 


1  1  +  2  log2 


and 


(  O 

W  /  _  \ 

1+sinJ 

1  (c0s2  )°8 

2 

6 

cos— 

2  J 

nin(n6)d6 


(nssl) 


was  evaluated  numerically. 

Equations  (All)  and  (A12)  constitute  N+l  linear  equations  for  A  and 
an  (n«l,2,3,...N).  The  formula  (12)  for  X  works  out  neatly  to 


-A-x^a. 


i-i 


■■1 


(A  13) 


and  the  independent  result  (10)  that  can  be  used  to  check  this  is  just 


X-U  -2A2/ar1/2 


(A14) 


i 
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An  alternative  procedure  that  gave  slightly  more  accurate  results  for 
a  given  number  of  simultaneous  equations  was  provided  by  the  use  of  the 
substitution 


H(s)  -  A(s/o)2  ♦  f(s) 


(A15) 


in  place  of  (A5).  The  first  term  in  (A6)  then  changes  to 


s 

2  + 


and  the  system  of  equations  (All)  and  (A12)  for  A  and  a.  (n«l,2,...N)  re¬ 
mains  valid,  except  that  Cm  and  Pq  change  to 


.  (3+a)Hm  a 


♦  =  Fm  +  =  W 


8 


m 


and 


Po.i 

“  10 


where 


H, 


4(m2-3) 

m(m2-4) 

4m 


m2-l 


(m  odd) 
(m  even) 


and 


Wm  *  (l-cos0)3log|-p^~jsin(m0)d0 


is  evaluated  numerically.  The  formula  (A13)  changes  to 


i-J 


■«  l 


(A16) 


$ 
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but  (A14)  stays  the  same.  For  values  of  a  up  to  30,  solutions  based  on  (A16) 
and  (A17)  with  N-60  were  judged  to  be  quite  accurate,  because  they  gave 
values  for  X  according  to  the  independent  formulas  (A  16)  and  (A14)  that 
were  in  excellent  agreement  with  each  other  (  1 %  discrepancy  at  a  >30,  and 
less  at  lower  a's).  The  final  results  shown  for  X  in  Table  I  were  obtained  by 
extrapolation  to  N>«>. 

Elastic-plastic  springs 

To  solve  (19)  for  g(s)  in  the  range  (0,ay),  we  differentiate,  as  before, 
and  again  use  the  substitution  (A3)  to  get 

OLy  (X  i— 

2sH'(s)-H(s)  +  sf  Mdt+gYf  —  -s  (A17) 

Jo  t-s  ^Jayt-s 

Now  with 

(A18) 
(A19) 

we  can  retain  the  Fourier  expansion  (A8)  for  f  (this  time  for  s  in  (0,  ay).  Once 
again  we  get  equations  (All)  and  (A12),  with  revised  formulas  for  Cm  and  Pq 
given  by 


and 


H(s)  -  A(s/oy)2  f(s) 


aY 

s  *  (1-COS0) 

2 


where 


The  expression  (21)  for  vP/vY  may  be  evaluated  as 


Vp/VY=-I  +  — . _^Og_=+l_J  >Og<.-»Y/a>  + 

-J 1  -^aY/a  +  y  (1  -a  Y/a)log  • 


l-JaY/a  2ab=1 


where 


rl+^aY/a  sin— 

cos  §■  sin(nS)  log - ■ - - 1-  d< 

-  2  sin| 


and  the  formula  (23)  gives 
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APPENDIX  B 

Consistency  of  Formulas  for  Toughening  Ratio 

We  want  to  show  that  if  g(s)  satisfies  the  integral  equation  (8),  then 
the  formulas  (10)  and  (12)  give  the  same  answer  for  X.  It  will  suffice  to 

prove  that 


N 


“s«aiT + 2g2 

'o  Vs  . 


+  2g2(a)  =  1 


if  g(0)  -  0,  and  Eq.(Al)  is  satisfied.  From  (Al) 


g(s)g'(s)  + 


ir \fim 

2  Jo  \  &  t- 


s  2  Vs 


and  integration  gives 


.ir*u 


To  show  that  (B3)  implies  (Bl),  we  just  have  to  verify  that 


fttg(s)ds 

Jo  Vs” 


f  2rr /T«dtds 

J  Jo  Jo  V  *  t_s 


But  if  the  integrals  exist  this  is  an  identity,  because 


■  Cf  ’SP**  •[(  £*f 


dsdt  =0 


l  «r.yvyvyv 
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APPENDIXC 
Elastic  Springs,  a  -»<*> 

If  we  let  x(a-s)/2  in  Eq.(Al),  and  then  let  o  ->  <*>,  we  get 


(UO) 


which  has  to  be  solved  subject  to  the  condition 


for  £  -4  <*> 


(Cl) 


(C2) 


that  is  consistent  with  the  far-field  stresses  (1).  Let  f(£)  denote  the  right 
hand  side  of  (Cl)  for  %  <  0,  and  say  g(£)  =  0  for  %  <  0.  Then  the  Fourier  trans¬ 
forms 

••  mm 

g (co)  *  J  g(S)e,0*d5  f(co)Hj  f(£)e,0*dS 

satisfy 


1^(i+io)i)g+=i[r+g(0)] 


(C3) 


Here  the  +  and  -  superscripts  denote  analyticity  in  the  upper  and  lower  half 
0D-planes,  respectively.  Now  split  (l+|ca|)  into  the  quotient  Q+/Q"  by  writing 

(Q*)  _  (CQ  _  <o/l(QI 
Q+  Q"  1+  |to| 


from  which  we  find 
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(Cfj  1 


-=L 


X/IXI 


Qf  2niJ-~  (1+|X|)(X-C0) 


■dx 


-  1  *  log  to 

2(1+®)  jt(co2-l) 


(Im(cD)  >0) 


Integration  gives 


Q+=  (const)V  1  +  co  expf--  f  iPi,zc*z 

l  Jo  z2-l 


Since  <d/|q>|  splits  into  (V <a)‘+/(vr  co)“ ,  Eq.(C3)  implies  that 

(Vco)+Q+g+  =  i[f  +  g(0)](V^)-Q- 


is  an  entire  function  of  to,  and  the  right  choice  will  turn  out  to  be  a  constant. 
Thus 


g(®)  = 


Vco(1+cd) 


exp 


i  f  log  zdz 

"Jo  z2_! 


(C4) 


the  only  singularity  of  which  is  a  branch  point  at  the  origin.  Evaluating  the 
inversion  integral 


g(^)  =  f  g(co)  exp(-io>5)  dto 
2  xw— 


via  integration  around  the  branch  cut  from  0  to  -i<*>  gives 


g(4)  =  Ce 


r 

,-«i/4  I 

Jo 


exp 


1  f  * log  T  dx 


1+x" 


0  VT(l+t2)3/4 


-dt 


(C5) 


For  large  £  this  can  be  evaluated  asymptotically  via  Watson's  lemma,  and 
matching  the  behavior  for  to  that  given  by  Eq.(C2)  gives 


[\y(3/2)-l  -log  £]  +  -•• 


gft) 


_1 _ 1 

(2k^)3/2 


( v[z]Er[z]/r[z]) 

with  C*  (l+i)/2.  To  get  an  expansion  of  g(^)  for  %  small,  use  (C4)  to  write  the 
Laplace  transform 


g(s)  m  I  e_s  *  g(5)  d£  =  g(is)  (s>0) 

JQ 


as 


g(s)  =  (2sf1/2(l+s2f1/4exp 


l_fs  log  p  dp 
***  1+p2  J 


Expanding  this  for  large  s  gives 


g(s) 


1  f  1  1 

vr[  s 


log  s  +1 

„2 


and  from  term-by-term  inversion  we  get 


g(4) 


08  9 +- 

y  =  -r(l)-  .57722  ... 


Finally,  with  Ce'*i/4  »  1  !•/  2,  Eq.(C5)  was  used  to  evaluate  g(£)  for  intermedi¬ 
ate  values  of  £  via  numerical  integration. 


TABLE  II  :  Experimental  data  and  material  properties 


